Mott Transition in Degenerate Hubbard Models: Application to Doped Fullerenes 
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The Mott-Hubbard transition is studied for a Hubbard model with orbital degeneracy N, using a 
diffusion Monte Carlo method. Based on general arguments, we conjecture that the Mott-Hubbard 
transition takes place for U /W ~ VN, where U is the Coulomb interaction and W is the band width. 
This is supported by exact diagonalization and Monte Carlo calculations. Realistic parameters for 
the doped fullerenes lead to the conclusions that stoichiometric A3C60 (A=K, Rb) are near the 
Mott-Hubbard transition, in a correlated metallic state. 
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The fullerenes have attracted much interest, not least 
because of the superconductivity of A3C60 (A=K, Rb) 
jl| . Much work has been performed to explain this phe- 
nomenon 0-0] starting from the assumption that the 
systems are metallic. The metalicity is an experimental 
fact ||, but surprising from a theoretical point of view. 
The Coulomb interaction U between two electrons on the 
same Ceo molecule may be a factor 1.5-2.5 larger than the 
width W of the (partly occupied) t\ u band [9|-[T5| . Since 
systems with U/W exceeding a critical value a bit larger 
than unity (~ 1.5) are believed to be Mott-Hubbard insu- 
lators [jlfj), one would expect this to apply to doped C60 
compounds. It has therefore been suggested that stoi- 
chiometric K 3 C 60 is a Mott insulator with a band gap of 
0.7 eV, and that experimental samples are metallic only 
because they are nonstoichiometric B . Alternatively, the 
critical value U c /W could be larger than commonly be- 
lieved. Which one of these alternatives apply, could have 
important implications for the properties of nonstoichio- 
metric A3C60 as well. We here study the Hubbard model, 
which has become the standard model for describing cor- 
related systems. In spite of the extensive studies of this 
model, U c /W is known accurately only for infinite di- 
mension in the absence of orbital degeneracy Jl6| . Since, 
however, almost all systems of physical interest have or- 
bital degeneracy |l7]] , we here focus on the effects of such 
a degeneracy, considering a three-dimensional fee lattice. 

First, we provide qualitative arguments suggesting 
that for a half-filled system with orbital degeneracy N, 
the critical value of U /W is enhanced by a factor y/N. 
These arguments are supported by exact diagonalization 
calculations for small systems and by diffusion Monte 
Carlo calculations for larger systems. The reason for the 
enhancement is that in the many-body treatment, the 
degeneracy increases the hopping rate of an extra elec- 
tron or hole in these systems relative to the one-particle 
band width. The degeneracy can thus stabilize a metal- 
lic state which may be highly correlated. This puts the 
doped fullerenes A„C6o (n =3, 4) (JV = 3) close to a 
Mott-Hubbard transition, so that a correlated metallic 
state of stoichiometric A3C60 is consistent with the large 



ond-mat/9608107 



values of U /W inferred from from experimental and the- 
oretical information. 

In a treatment using the Gutzwillcr Ansatz and the 
Gutzwiller approximation [fl8| for the Hubbard model, 
Lu found the critical value U/W to vary as (N + 1) at 
half- filling p9| . Here we provide qualitative arguments 
showing that the degeneracy is important, but that the 
critical value is substantially smaller and has a weaker 
(VN) dependence than in the work of Lu. Our argu- 
ments are supported by exact diagonalization and Quan- 
tum Monte-Carlo calculations. 

To discuss the Mott-Hubbard gap for a model with the 
orbital degeneracy N, we first consider a (bipartite) lat- 
tice with M sites and a hopping integral t im j m > = tS mm ' 
between orbitals with the same quantum numbers m on 
neighboring sites i and j. For half- filling the energy gap 
is then 

Eg = E(NM + 1) + E(NM - 1) - 2E(NM), (1) 

where E(L) is the ground-state energy for L electrons. 
For £/>twe have 

E(NM) = ^N(N-l)MU + 0( t -), (2) 

since no electron can hop without an extra energy cost 
U . For instance, this energy is obtained for a Neel state 
|0), with the moment N. In the system with NM + 1 
electrons, the extra electron can hop without an extra 
energy cost U. We thus consider the states 



|1>=<J0) 

10 = -T^iuX^imT^Tl )' 



(3) 



where V'Jmj creates a spin up electron on site j in orbital 
to. The sites i are nearest neighbors of site 1, which we 
have assumed to have spin down electrons in the state 
|0). We then have 



(i\H\l) = VNt, 



(4) 



i.e., an electron hops from site 1 to site i with a matrix 
element which is a factor \/~N larger than in the one- 
particle case [pOl . This is due to the fact that any of the 
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N spin up electrons can hop from site 1 to i. We next 
construct states \ij) (j ^ 1) where a nearest neighbor j 
of i has an extra electron. The corresponding matrix el- 
ement is y/~Nt. However, the states \ij) depend on which 
site i had an extra electron in the intermediate state. We 
use the analogy with the one-electron case, although it is 
not rigorous, since in the many-body case there is a string 
of reduced spins along the path of the moving additional 
occupancy. In the one-electron case, the hopping of an 
electron lowers the energy by £5, where £& is the bottom 
of the band. The close analogy between the one- and 
many-particle problems then suggests 

E(NM + l)xiE(NM)+NU + VN£b, (5) 

with an extra factor y/N in front of Eb- Using a similar 
result for E(NM — 1) we obtain 



E g = U- VN{e t -s b ) = U- VNW, 



(6) 



where £t (= — £& in this model) is the top of the one- 
particle band. This can be generalized to the case where 



all the hopping integrals t 



— t for i and 



j near- 



est neighbor sites. For N = 1 we obtain E g = U — W. 
If this result is extrapolated to smaller U, it predicts a 
Mott-Hubbard transition for U/W ~ 1, close to previous 
results. Eq. (g) suggests, however, that for a degenerate 
system the transition takes place for a larger ratio ~ \f~N . 

To test this, we have performed exact diagonalization 
calculations for a Hubbard model with a six site (3 * 2) 
lattice with (N = 2) and without (N = 1) orbital de- 
generacy. Table fl compares the exact band gap with the 
estimate U — y/NW . The exact results agree very well 
with the simple estimate for large U, confirming the es- 
timate in the limit it was made. For intermediate U the 
agreement is less good. The reason could be that the gap 
is finite for a finite system or that the extrapolation of 
the estimate to intermediate U is invalid. To check this, 
we consider larger systems. 

To study A3C60 we have considered a model, which 
includes the three-fold degenerate, partly occupied t\ u 
level. We include the on-site Coulomb interaction U 
and the hopping integrals t im j m > between the molecules. 
This leads to the Hubbard-likc model 



H = E E ^ 

ia in— 1 

+^E E 



(7) 
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where the sum < ij > is over nearest neighbor sites. 
We have used a fee lattice. The hopping integrals U m j m i 
have been obtained from a tight-binding parametrization 
pT| , p2| . The molecules are allowed to randomly take one 
of two orientations in accordance with experiment p3[ , 
and the hopping integrals are chosen so that this ori- 
entational disorder is included fl241. The band width of 



the infinite system is W = 0.63 eV. The model neglects 
multiplet effects. The inclusion of these effects may fa- 
vor antiferromagnetism, which could reduce the critical 
value U c /W for a Mott-Hubbard transition. 

We first consider a cluster of four Ceo molecules, for 
which the Hamiltonian (^) can be diagonalized exactly. 
The results are shown in Table | for the cases of three 
(A3C60) electrons per site. Except for small values of U, 
the gap Eg is actually smaller than than U—y/SW, which 
may be due to the hopping being more complicated than 
discussed above. The importance of the degeneracy is, 
however, clear. 

To study the Mott-Hubbard transition, we consider 
larger systems. We use a diffusion Monte Carlo method, 
which we have developed along the lines of ten Haaf et 
al. (25). Starting from a trial function \SS?t), we calculate 



= [i-t(h 



l |* T > = F n \* T ), 



(8) 



TABLE I. Band gap as a function of U and iV compared 
with U - y/NW for a six-site model (N = 1 and N = 2) 
and a cluster of four Ceo molecules. W=0.6 for the six atom 
cluster and 0.58 for the Ceo cluster. E g (n) is the band gap 
for n electrons per site. All energies in eV. 

U N = 1 (6-site) N = 2(6-site) N = 3(4C 60 ) 

E g {l) U-W E g {2) U-sf2W E g (3) U - s/%W 



10 
5 
3 
2 
1 



9.35 
4.39 
2.44 
1.48 
0.54 



9.40 
4.40 
2.40 
1.40 
0.40 



9.06 
4.12 
2.18 
1.24 
0.35 



9.15 
4.15 
2.15 
1.15 
0.15 



8.39 
3.48 
1.64 
0.90 
0.46 



9.00 
4.00 
2.00 
1.00 
0.00 



TABLE II. The total energy for four Ceo molecules with M 
electrons and U = 1 according to the variational (VMC) and 
diffusion (DMC) Monte Carlo calculations compared with the 
exact results. The Gutzwiller parameter is given by g. The 
statistical error is estimated to be two units in the last digit. 



M 


g 


VMC 


DMC 


Exact 


11 


0.45 


8.4970 


8.4677 




11 


0.50 


8.4842 


8.4677 


8.4649 


11 


0.55 


8.4873 


8.4684 




12 


0.45 


10.7364 


10.7014 




12 


0.50 


10.7226 


10.7009 


10.6994 


12 


0.55 


10.7249 


10.7012 




13 


0.45 


13.4363 


13.3990 




13 


0.50 


13.4211 


13.3988 


13.3973 


13 


0.55 


13.4233 


13.3991 




Eg 


0.45 


0.4605 


0.4639 




E,j 


0.50 


0.4601 


0.4647 


0.4634 


Eg 


0.51 


0.4608 


0.4651 
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where w is an estimate of the ground-state energy. If r 
is sufficiently small and \^t) is not orthogonal to the 
ground-state, |<J/( n )) converges to the ground-state as 
n — > oo. |\I>( n )) is calculated in a Monte Carlo approach. 
Although this approach in principle is exact, it suffers 
from the usual sign problem. We therefore use the lat- 
tice equivalent of a fixed node approximation pjJ, and 
introduce an effective Hamiltonian for which the nondi- 
agonal terms are zero if 



(^ T \R) (R\f\R'} (R'I^t) < 0, 



(9) 



where R represents the coordinates of all the electrons. A 
term is added to the diagonal part of the effective Hamil- 
tonian, so that its energy is an upper bound to the exact 
energy, with the upper bound agreeing with the exact 
energy if is the exact wave function [ p5| . 

For the trial function we make the Gutzwillcr Ansatz 



I*t> = s»*|*o> 



(10) 



where g is a Gutzwiller parameter and the power 
n,j is the number of double occupations rid = 
X)i Y,am«T<m> n i<ymn ia ' m ' in the Slater determinant |* )- 
To construct |^o) we solve the Hamiltonian (]?]) in the 
Hartree-Fock approximation, replacing U by a variational 
parameter Uq- Depending on Uq the trial function is 
paramagnetic or antiferromagnetic. 

As a test, we study a cluster of four Ceo molecules. 
We calculate the expectation value of H for \^t) in a 
variational Monte Carlo (VMC) approach and perform 
a diffusion calculation (DMC) . The results are compared 
with the exact results in Table ||. The agreement is very 
good, with an error of 0.003 eV or less in the DMC cal- 
culation for the parameters considered. 

To obtain an efficient extrapolation to infinite systems, 
we correct for finite size effects . For the present prob- 
lem, there should be a gap for a finite system even for a 
small U . We can see this by distributing the charge uni- 
formly over all the levels and calculate the electrostatic 
energy. This gives a contribution U /M to the gap |27| , 
where M is the number of molecules. We then introduce 



E g = E g -^-E g (U = 0), 



(11) 



where E g (U = 0) is the gap in the one-particle spectrum. 
We may then expect E g to be rather independent of M 
for itinerant systems. Both U/M and E g (U = 0) go to 
zero for M — » oo, so the subtracted quantities improve 
the convergence but do not change the M — ► oo result. 

In Fig. [l] we show results for E g that suggest a Mott- 
Hubbard transition for U between 1.5 and 1.75 eV. This 
corresponds to a critical ratio U c /W ~ 2.5. This is sub- 
stantially smaller than the ratio 4 found by Lu M| • F° r a 
fully frustrated system with infinite dimension but with- 
out orbital degeneracy U c /W — 1.5 was obtained p6[. 



Multiplying by the degeneracy factor y/3 leads to the ra- 
tio \/3 x 1.5 = 2.6 in good agreement with our results. 
It is interesting to note that for U=1.75, which gives an 
insulator, an antiferromagnetic trial function results in 
the lowest energy for the larger systems. 
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FIG. 1. Eg (Eq. ( |ll[ ) for different U and as a function of 
1/M, where M is the number of molecules. The band width 
varies between 0.58 eV (M = 4) and 0.63 eV [M = oo). 

The inclusion of higher subbands may increase U c /W 
further. In the large U, half-filled case, the hopping to a 
higher subband costs an extra "diagonal" energy U + As, 
where Ae is the energy separation of the bands. The cor- 
responding energy cost in the states with one extra elec- 
tron or hole is only As. The hopping to higher subbands 
should therefore lower the energy more in the latter two 
states, which decreases the gap. Simultaneously, the one- 
particle width of the t\ u band due to hopping into higher 
levels on neighboring molecules is reduced by a few per 
cent p2|| . Only bands with Ae < U contribute to the 
reduction of the band gap. A simple model suggests that 
this could increase U c /W by 0.1-0.2 for A^Cqq. For 3d, 
4/ or 5/ compounds (TV = 5 or 7), a crystal-field split- 
ting may lower the degeneracy. The higher crystal-field 
split states cannot, however, be neglected when deter- 
mining U c /W, unless the crystal-field splitting is much 
larger than the Coulomb interaction. 

The average Coulomb U has been estimated to be 1.6 
eV H and 1.4 ± 0.2 eV jl(| from Auger experiments for 
the occupied orbitals. A slightly smaller (0.2 eV) value 
was found for the highest occupied orbital (h u ) M, and 
U may be somewhat further reduced for the higher-lying 
ti u orbital. The reduced screening in the surface layer 
is estimated to increase U in this layer by 0.3 eV [|. 
The surface sensitivity of Auger may therefore slightly 
overestimate U. Theoretically, U has been estimated to 
be in the range 0.8-1.3 eV [0. Band structure calcu- 
lations have found the band width 0.41 eV Jl3]] , 0.6 eV 
H and 0.65 eV (m) for K 3 C 60 . Many-body effects were 
found to increase the band width by about 30% in the 
GW approximation for undoped Cgo |15| ]. The many- 
body effects treated in the Hubbard model, due to the 
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interactions in the t\ u band, are not included in the GW 
calculation, since the t\ u band was empty in the GW cal- 
culation. Using the GW enhancement of the band width 
is therefore not a double-counting. For RbsCgo the band 
width may be reduced by about 15% p2[ , compared with 
K 3 C6o- This leads to a likely ratio U/W in the range 
1.5-2.5. In view of our QMC results, this is probably too 
small to cause a Mott- Hubbard transition for A 3 Ceo- 

One could ask what are the effects of U > W in a 
metal. The system can remain a metal in spite of the 
large U/W thanks to the hopping via several channels 
(Eq. (H)). As in the impurity problem |2(J, this may fa- 
vor the formation of a singlet state. In A3C60, which 
has an odd number of electrons per site, the singlet is 
consistent with a correlated non-magnetic state, which 
must be metallic if it obeys the Luttinger theorem p8[ . 
It is suggestive that such singlet correlations could tend 
to promote superconductivity in these systems. 

While A3C60 are metals, A4C60 are insulators p9| , 
although band structure calculations predict that they 
should be metals |$0|. A4C60 is therefore not a band in- 
sulator, and this raises interesting questions about the 
difference between A3C60 and A4C60- The presence of 
four electrons per molecule in A4C60 favors the forma- 
tion of a static Jahn- Teller effect, which may open up a 
gap |3l) , in particular if electron correlation effects are 
also included. The electron-phonon coupling should re- 
duce the effective hopping for both A3C60 and A4Cgo- 
One may also speculate that the gap in A4C60 is a conse- 
quence of the interactions in a case with an even number 
of electrons per site. This could be a correlation induced 
gap with no development of magnetic moments, i.e., a 
different type of order than in other correlated Mott in- 
sulators such as NiO. Besides A4C60 is bipartite while 
A3C60 is not, which may favor the development of or- 
dered phases (anti-ferromagnetic or other types of order) 
and the opening of a gap in A4C60- Finally, the Fermi 
surface of A4C60 has almost perfect nesting, which may 
favor a charge density wave and the opening of a gap |j0| . 

To summarize, we have shown that the orbital degener- 
acy increases the value of U/W where the Mott-Hubbard 
transition takes place, due to the more efficient hopping. 
This puts the critical ratio for A3C60 at the upper range 
of ratios U/W estimated for doped Ceo compounds and 
makes it likely that A3C60 (A=K, Rb) is on the metallic 
side. Although experimental samples of A3C60 may be 
nonstoichiometric, it is therefore not necessary to assume 
so to explain the metalicity. The large value of U/W may 
nevertheless lead to substantial correlation effects, which 
may have a large influence on certain properties. Most 
systems of interest have orbital degeneracy |l7t| . There- 
fore our results should be relevant also for many other 
systems close to a Mott-Hubbard transition. 
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